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\Q • A dynamical symmetry for spherical collapse has been studied using a linear transformation of 

the initial data set (mass and kinetic energy function) and the area radius. With proper choice of 

, the initial area radius, the evolution as well as the physical parameters namely energy density and 

shear remain invariant both initially and at any time instant. Finally, it is found that the final 

' outcome of collapse depends on the initial choice of the area radius. 

Debnath, Chakraborty and Dadhich in a recent paper [1] while demonstrating the true invariance of the 
collapse dynamics, chose to scale the comoving radial coordinate also along with other initial data sets. This 
has, in fact, destroyed the invariance of energy density and shear as well as the collapse dynamics itself. It can 
be shown clearly by the following arguments. In their case the linear transformations are chosen in the form 



> 

(N 



X 



PACS numbers: 04.20 Dw, 04.20 Ex 



F -> a"- L F, / -> cr/, R->a R and r -> a r, 

where a is a constant, n is the total number of dimensions and r is the comoving radial coordinate. F, /, R 
\& have the usual meanings [2,3]. The new radial coordinate in reality changes the magnitudes of the functions 
■ F, f and R. So when the above scalings are performed the physical quantities such as p, a etc although 
""^ \ retain their structural forms unchanged in terms of F, /, r etc they do suffer changes in their magnitudes 
O"^' an d can no longer be claimed invariant. The scaling of the coordinate r was, however, demanded by their 
initial choice R(U,r) = r. One can avoid such unwanted scaling of r if one chooses R{ti,r) as another 
bX). suitable function of r. In the following text we have done exactly the same. So our choice gives the true 
invariance of p, a and the dynamical equations for a class of data sets which lead to the same end results as 
well as the evolution of the system chosen. We should point out here that though our model is a spherical 
Tolman-Bondi-Lemaitre (TBL) model yet similar results may also be valid in the quasi-spherical collapse as well. 



The inhomogeneous spherically symmetric dust space-time of n dimension can be represented by the following 
metric [4] 

ds 2 = dt 2 Trrdr 2 - R 2 dSl 2 n 2 (1) 

1 + /0) 

Here we have used comoving coordinates. In the above metric goo (the co-efficients of dt 2 ) = l, because the 
dust particle follows geodesies. The function f(r) is related with the curvature of the (n— 1) dimensional space 
part and fixes it as bound, unbound and marginally bound for f(r) < 0, f(r) > and f(r) — respectively. 
The function R(t, r) stands for the area radius of the corresponding shell. The evolution equation is given by 



R 2 = fir) + |S (2) 
where F[r) can be interpreted as the mass contained within the comoving radius V and is reflected in the 
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expression for energy density p(t, r): 



, . (n-2) F' 



The shear scalar cr(t,r) is given by [3,4] 



n-2 {RR'-RR') 

a= \j2(n—T) RR> (4) 

The above collapsing models are characterized by the initial data given on the initial hypersurface t — U. It 
is to be noted that one can make an arbitrary relabeling of spherical dust shells by r — > ip(r), without any loss 
of generality. We make here the suitable choice of the radial coordinate such that at t = ti, 

R{U,r) = — (5) 

where, a and b are constants. 

We now proceed to find out a class of initial data functions which under linear transformation yield exactly 
the same density, shear and evolution equation not only at the initial instant but also at any arbitrary instant. 
These are [1,5] 

{F, /} -> {mF, If} (6) 

In order that evolution equation (2) remains invariant under the above transformations we must require a 
relation between m and I in the form 

to = (7) 

and a scaling of the area radius as 

R^l l ' 2 R (8) 

So finally we find that the transformations (6) and (8) in view of (7) retain p and a invariant. The interesting 
point is that the above transformations do not alter the scaling used for the radial coordinate given in (5) 
provided the constants a and b are related as follows 

a(n-l)-26=l (9) 

It is now expected that the collapse of the dust sphere would start from a regular initial hypersurface t = ti 
and hence pi and cr, should be regular on the initial hypersurface and in particular at r = 0. 



The initial density and shear can now be given using (5). These are 

p(n-l)o-l r pf ft 



. . (n-2) rF' 
Pi = p{U,r) = 



2 F 
and 



r^^^(l + ar^-br f) (i + r »-8^^) 

(11) 

Now if we choose a power series function for F(r) and f(r) around r = 0, we may conclude that both 
and ^j- arc finite at the centre and hence in order that p(t i; r) is regular at r — > 0, we must have 



Urn r n 1 F< w = a non-zero finite number. 
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The regularity of the initial shear at r — > further demands in view of (11) that the quantity 

r (n-3)o-l 



Urn r n 3 — 
r -> 



/( „-3) 6 



to be zero or a finite constant. 



Since our study is restricted to the region close to r = for smooth initial data we choose the power series 
functions for F(r) and f(r) in this region, so that we can write [4] 



3=0 



(12) 



(13) 



J=0 



Using the expressions (12) and (13) close to the centre r = 0, the condition worked out previously for the 
regularity of the initial density demands a and (3 to satisfy the following relation 



aa — b/3 + 1 



a 



n- 1 



(14) 



We now proceed to examine the expression for shear given in (11). In view of (12) and (13) and also (14) 
the quantity -p — (n — 1) ^1 + ar ^ — br reduces to zero. So in order that cr(tj,r) vanishes at the centre 



(r = 0) 



-3 F { 



0. 



f( „- 3 )b-i must vanish at r 
The above in turn leads to the condition 

(n - 3)(1 + aa - b/3) > a - (i, 

which when combined with the condition (14) imposes the following restriction on the magnitudes of a and 
(3 namely, 



> 



a- (3 



(15) 



provided we assume a > (3 and n > 3. 

So finally, we get the restriction in the form 



a > [3 > 



2a 



(16) 



We now proceed to write the expression for the initial energy density near r = using the power scries 
expansion forms (12) and (13). It is given by 



Pi = p{U,r) 



6(71-1) 



(n-l)(rc-2)/° 

2 F 2fc 



(17) 



where p2 is actually a constant, the exact expression of which is given later. 



The equation (17) shows that for a smooth physically acceptable behaviour of the density function at r = 



one can demand ( ) = and hence it follows that the constant 

V / r=0 

JO ^0 

The constant pi in (17) has the simplified expression using (18) in the form 



Pi = b 1 + 



(n-1) 



^-^1+26^' 



1_ {2-2a + 6(n- l)a} 
a (n — l)a 



(18) 



(19) 
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Again if one assumes that pi is maximum at r — and decreases with increasing r one must have p- < in 
the region close to r = 0, which in turn demand p2 < near the centre. 



Now the solution of the evolution eq(2) can be written as [4] 
2 



t-U = 



{n-l)VF 



' n^F^ r l F a(„-3)-l 



n-1 1 fR n ~ 3 

T ]-R— 2 F 1 [-,d,d + l,- J 



(20) 



where d = i + and 2 fi is the usual hypergeometric function. 

If t = t s (r) stands for the time of collapse of a shell of radius R giving rise to a shell focusing singularity at r, 



2 n-l F ai " ^ 1 EV»(fl-3)-l 

*.W - = — x r ^ jh^h 2i % d > d + l ' ~ r /6 (n- 3) -i 1 



So the time of central singularity is given by 



2F b 



(21) 



(22) 



t = Urn t s {r) =U + - , 

r^O (n-l)/ 2 

Here the power series expansions (12) and (13) have been used. 

If t = t a h{r) stands for the time of formation of apparent horizon at coordinate distance V then area radius 
is restricted by [3,4] 



and we have from (20) 



tah(r) - U 



n-l 



r 2 



«("-D-i 

F 2 



R(t ah (r),r) = [F(r)} — 

pa(n-3)-l 2F^ 



(23) 



.1 



„n— 3 J 



j, b(n-i) 2F 1 [-,d,d+l,-r /fc( „_ 3) _ 1 



]- — ^[-,^+1,-/] (24) 



Therefore, using the series expansions (12) and (13) the time difference between central shell focusing singu- 
larity and the apparent horizon at coordinate distance r, (close to r = 0) is given by 



tah(r) - t Q = 



F~ 



(3n-7)q-3 



(n-l)(d + l) 
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r 2 + 0(r 3 ) 



(25) 



A little analysis of the equation (25) shows that when b > 1/2 we have 1/6 < 2, so that there is only a black 
hole and no naked singularity. The magnitude of is 3 at n = 4 and decreases with the increasing number 
of dimensions. For b < 1/2 and n < 5 there may occur either a black hole or a naked singularity depending on 
the sign of the coefficient of r 2 in equation (25). Again for the number of dimensions greater than five, i.e., for 
six dimensions or more, < 2, so that we get only a black hole. This conclusion is exactly identical with 
that obtained in our previous paper [2] . 



by 



Marginally bound case (/ = 0): 

Here our initial choice is R = (r + r )F a , r and a being positive constants. The density and shear are given 



, , (n-2)F'(r) , , / (n - 2) 



R' 



R 
R 



(26) 
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and the evolution equation simplifies to 



R 2 = 



Fjr) 



(27) 



Suppose we perform linear transformations on F and the area radius R in the following manner: 

F->pF, R->qR (28) 

where p, q are constants. The invariance of p, a and the evolution equation (both initially and at any time 
instant) under the above transformations yield 



n-l 



and p = q 



n-l 



On the initial hypersurface, the energy density and the shear scalar are given by 

(n - 2) 



rF' 
F 



and 



2 (n-2) 

CT; = ~ 



(n - l) 3 r 2 



(r + r ) n+1 



{n-Vjr . r F 



r+r 



+ — 



(29) 



(30) 



(31) 



Here it is clear from the above that ro ^ 0, because otherwise there will be a permanent singularity of infinite 
energy density and shear at the centre. In this case the solution of the evolution equation has the simple form 



t-U = 



(n-l) 



(r + r )" 



(32) 



So proceeding as before 



t a h(r) -t =r Q 2 r 



n-l 



Fo~ 3 r— , 



(33) 



which shows that t a h{r) > to for small r (i.e., close to singularity). Hence we always have naked singularity in 
any dimension. Finally, it is to be noted that the parameter ro plays a crucial role in characterizing the end 
state of collapse. 



In conclusion we must mention that a significant point in our paper is that we have avoided choosing 
R(ti,r) — r as is done in all the previous papers without perhaps any exception. Since R(ti,r) can be chosen 
any arbitrary function of r without any loss of generality we have preferred to select a different function of r 
for R(ti,r), which allows us to perform a different set of scalings and consequently draw different conclusions 
in some situations. For example, the sign of the quantity (t a h{r) — t ) in our present work depends upon an 
additional parameter b, the suitable choice of which may yield in certain cases different conclusions regarding 
the occurrence of the naked singularity or black hole. However, the result we obtained in one of our papers 
that there is no naked singularity in the space-time having more than five dimensions, still holds showing the 
possibility of this result being valid in general. 

It is to be noted that in our paper even if the quantity (t a h{r) — to) is not invariant, its sign, however, will 
remain unchanged under the scaling chosen in our paper, as one can easily verify from the equation (25). So 
the conclusion about the end result of the collapse (naked singularity or black hole) still remains valid even 
after the scaling the parameters. Since it depends on the sign and not the magnitude of the quantity (t a h(r) — to). 

It is further to be stresses that for the parameter b under certain restriction (b > 1/2) the end state of the 
collapse always leads to black hole irrespective of any number of dimensions of our space-time. So finally we 
may conclude that with a suitable choice of the initial area radius R as a function of the comoving coordinate 



G 



r we can ensure the validity of the so-called Cosmic Censorship Conjecture. 
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